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We consider Andreev reflection in ferromagnet-superconductor junctions and extend the treat- 
ment by de Jong and Beenakker to incorporate the effect of a repulsive potential at the inter- 
face. We show the existence of a spin accumulation at a finite voltage in this model. 



1 Introduction 

The interplay between spin polarization and superconductivity has recently gained a renewed 
interest because of the nowadays possibility to study Andreev renectipn0 in small junctions in= 
volving ferromagnetic metals in ballisticBH as well as diffusive systemsQ. de Jong and BeenakkerH 
studied multichannel effects in Andreev reflection in ferromagnet-superconductor junctions with 
a transparent interface. Fal'ko at ai.Q and Jemeda et aJ.0 have recently shown the existence of 
spin accumulation in a diffusive ferromagnet-superconductor interface model. We consider here 
ballistic ferromagnet-superconductor interfaces and extend the de Jong and Beenakker treat- 
ment B to incorporate a finite repulsive interfacial potential. We show the existence of a spin 
accumulation also in this model. 

2 One channel FS junction model 

We first consider a one channel model describing a contact between a ferromagnet and a su- 
perconductor. Both are thermally equilibrated reservoirs. Their energy distribution is therefore 
given by the equilibrium Fermi-Dirac distribution. We consider a Stoner ferromagnet, modeled 
by non interacting electrons in the presence of an exchange field h(x) = h9(—x). The polarization 
enters this model via the ratio hjep of the exchange field h over the Fermi energy €f. We also 
consider the existence of an interfacial repulsive Dotential H5{x) at the interface. We introduce 
the dimensionless barrier strength Z = H/hvpU. For a metallic junction Z = and Z = +oo 
for a tunnel junction. The superconductor is a BCS superconductor with a gap A(x) = A9(x). 
The scattering problem will be solved following Blonder, Thinkham and Klapwick (BTK)Q. We 
consider an electron with an energy E incoming on the junction. Since this electron also has a 
spin energy term, its wave vector depends on spin polarization: 

Hk^ = V^rn^E +Jb±h. 

Whereas de Jong and Beenakker considered multichannel effects with Z = oi, we focus here 
on the case of a finite value of Z, similarly to Zutic and Vails 13. Let us consider an incoming 
electron with a spin up and match the wave functions on the left and the right of the F-S contact 
at a constant energy E. The left-side wave function is 

$ L {x) = ^ J e lk r x + a e \z\Z l ,E) ^ e lk r x + b e \z\ Z l , E) ( J e~ ik > 

while the right-side wave function is 

1> R {x)=c( U0 )e ik ^ + d( Vo )e~ ik ^, 



where the superscript "e j" refers to an incoming electron with a spin up and no and vq denote the 
usual BCS coherence factors. The Andreev reflection amplitude a, the backscattering amplitude 
b and the transmission amplitudes c and d are found by matching the wave function and its 
derivative according to 



MV = Mo) (i) 



This leads to 



l + 6 e T = CUo + duo = ^( CUo _ duo )__l_( 1 _ & eT) ( 3 ) 



i e] ' = cv + ddu = ^(cv - du ) - ^i ^ , ( 4 ) 

with a different barrier Z 1 ' = mH/h 2 kp and = mH/h 2 kp for spin up and spin down electrons. 
Solving for the Andreev reflection coefficient, we find 

cf\z\zKE) = 7 — , 4U0V0 — (5) 

{(«8"«g)(5 + 4ZZT + ^)+l + |t + 2i( U §-t,g)z(|t-l)} 

This expression coincides with the BTK Andreev amplitude! if Z^ = Z^ = Z. 

Let us now consider an incoming hole in the spin down band. This will be Andreev into 
an electron with a spin up, therefore contributing to the current of left-moving electrons with a 
spin up. The left-side wave function is 

$ L (x) = ^ J ^ e~ ik r x + a hl (Z\Z l ,E) ^ J ^ e~ ik ^ x + b h \z\ Z l , E) ^ J ^ e lk r x 



and the right-side wave function is 



^ R (x) =c^ V u ° o J e~ lkFX + d y v ° J eikFX - 

The corresponding matching equations are 

a h \z\z\E) = cv + du = -^(cv -du Q ) + ^- f a h t(Z\Z±,E) (6) 

l + b hl (ZlzKE) = cuo + dv = -^(cu -dvo) + ^ l (l-b h ^Z\zKE)). (7) 

Comparing Eqs. |6] and ^ with Eqs. |3| and |], the Andreev coefficient a h ^ (Z* ,Z\E) of an incoming 
hole with a spin down is deduced from Eq. [5] with the prescription Z — ► — Z, Z^ — ► — Z^ 
and ZJ- — ► — Z*. Since the incoming and outcoming wave vectors have been interchanged, we 
interchange Z^ and Z\ and we change sign because we are using hole instead of particle wave 
functions. We obtain a h ^(Z\ Z^,E) = (k^ /k^)a e ^ (Z\ Z^,E). This identity is expected because 
of particle-hole conjugation. 

Now an incoming electron with a spin down has an Andreev reflection amplitude 

a el {Z\Z l ,E) = a e \Z l ,Z\E), 

which is not related in a simple way to a e ^ {Z\ Z^ , E) below the superconducting gap because 
of the phase shift term in Eq. |5|. A final set of identities is given by the amplitude current 
conservation, already discussed by Zutic and Vails 

fc T = k L \a e1; {Z\Z l ,E)\ 2 + k)\b e \z\z\E)\ 2 . (8) 



Let us now discuss the consequences for the zero temperature current below the superconducting 
gap. We find 

/TOO = /*M*-.V)-/*$)W,ZW*-M-Z-.V)) W 

- j dE\b e \zlzl,E)\ 2 fo(E-eV) (10) 

= JdE(^\a^(Z\z\E)\\UE-eV)-f (E + eV)), (11) 

where /q = 9(—E) denotes the zero temperature Fermi distribution. Similarly, 

I l (V) = JdE (j^j \a^{Z\ ZKE)\\f {E - eV) - f (E + eV)). 

The current of spin down electrons is not equal to the current of spin up electrons because 
of the phase shift term in the expression Eq. || of the Andreev reflection coefficient, there- 
fore leading to a spin accumulation at finite voltages. If we consider the limit of a trans- 
parent interface (Z = 0, Z^ = and Z^ = while Z^ jZ and Z^ jZ remain finite) we find 
a e ^-{Z\ Z^, E) = (k+ /A;T) a e t (Z\ Z^ ,E), showing that the current is not spin polarized if the 
interface is transparent. 



3 Spin accumulation 



Since a spin polarized current cannot enter below the superconducting gap, we should consider 
a mechanism restoring a spin unpolarized current at finite voltages. We consider spin flip 
processes in the ferromagnet giving rise to a different chemical potential for spin up and spin 
down carriers at the F-S interface: <fyJ = —5^ = 5[i. We will solve for 5a and show that 
8fx _> 0. This existence of a spin accumulation was put forward by Fal'ko et al.u and Jedema et 
al.u in the diffusive regime and we show the existence of a spin accumulation also in the ballistic 
regime. The currents of spin up and spin down electrons are 




{Z\ Z\ E)\ 2 [f (E -eV- Sfi) - f Q (E + eV + 6»)} (12) 



(Z\ Z\ E)\ z [f (E -eV + 5fx) - f (E + eV - Sfi)] , (13) 



Physically, an outcoming electron with a spin up can be produced by backscattering or Andreev 
reflection. Backscattering is controlled by the chemical potential +S/J, of the spin up electrons 
whereas Andreev reflection is controlled by the chemical potential —6/j, of the spin down electrons. 
The chemical potential 8jj, is then determined by imposing a zero spin current entering below the 
superconducting gap: E(6fj,, V) = I+(Sfj,,V). We solved this equation in the small-//, small-y 
regime, in which case the integrals in Eqs. [l2| and [D] can be expanded in the small parameter 
eV + 5/j,. We note (E) = (k+ /k^)\a e ^ (Z\ Z\ E)\ 2 and make the following approximation for 
the integral Eq. Ol of the current: 



g \E)dE^ l -(g\0)+g\eV + 5^) (eV + Sf,). 



We obtain 



7 l 



5 T (0) +g\eV + £//)] (eV + Sn) 
g l (0) + g l (eV - 6fi)] (eV - Sfi). 



(14) 

(15) 
(16) 



Z=1 




Polarization 

Figure 1: Variation of the chemical potential <5/i = Sfi^ (0) = —6^(0) as a function of spin polarization for Z = 1. 
The chemical potential is measured in units of the superconducting gap A = 1. The voltage is V = 0.1. 



If only the first terms 5^(0) = g^(0) had been included, one would have 5fi = whereas we find 
a finite and positive 5/j, when the second terms are included. If V is finite and 5fi = 0, a spin 
current is present (jf < 1^) and the value of S/i is determined by equating the spin up and spin 
down currents. The variations of <5/i as a function of the spin polarization P are shown on Fig. |l|. 



We have approximated the integral Eq. |12] according to Eq. 14. This approximation does 
not modify the qualitative behavior because g^(E) and g^(E) are increasing with energy. Even 
without this approximation, the spin up current increases upon increasing Sfi, and the spin down 
current decreases, therefore leading to a finite 5fi at finite voltages. 

A finite value of 5(i at the interface can be generated by spin flip processes in the ferromagnet. 
If we note f R (E, x) = fx(E— ji° (x)) the distribution of right moving electrons in the ferromagnet 
(x < 0) the semi classical Boltzmann equation with a spin flip collision term is 

§~JU E ^) = -Y s {fl(E,x)-f l R (E,x)) (17) 
^fk(E,x) = ~(f l R (E,x)-f R (E,x)), (18) 

with l s the spin flip length and the boundary condition 5fj,'(0) = —5^(0) = —5fi. From what 
we deduce dfj, ± (x) = ^e 2x ^ a 5fi. 

4 Multichannel model 

Let us now consider multichannel effects discussed by de Jong and Beenakker in the case of a 
metallic contact B. We consider a number iV-'- ~ (k^-) 2 of spin down channels, and a number 

iV^ ~ (k^) 2 > A^-^ of spin up channels. We consider all the spin up and spin down channels to 

t i 

have the same Fermi wave vector k F and k F . The transport equations below the superconducting 
gap are f R (E) = N° f (E - eV), and 



fl(E) = NUl-r) \a h ^{Z\ Z^E)\ 2 [l - f (-E - eV)} + Z\E)\ 2 h{E - eV) 




+(iVT - N l )f (E - eV) 

,„.,((!;),..„«.„,,-,,-.-.„,*,,.«,.,.<.-«). 

the last term arising from the fraction 1 — N^/N^ of the incoming electrons with a spin up that 
are not Andreev reflected and are fully backscattered. The zero temperature current of spin up 
and spin down electrons is 

„ = Hif^y^zWWB-^-Ml + 'V)) (») 

I 1 = N l JdE(^j\a ei (Z\Z^E)\ 2 (ME-eV)-f (E + eV)). (20) 

This expression is equal to the single channel current times the number of spin down channels. 
A similar mechanism of spin accumulation also exists in this model, the spin accumulation being 
independent on the number of channels. However 5p is small even though finite and therefore 
should not give a significant contribution to the junction resistance. 

5 Conclusions 

We have analyzed transport in resistive F-S junctions and shown the existence of a small, even 
though finite spin accumulation at finite voltages. A spin accumulation was obtained by Fal'ko 
et aJ.Q and Jedema et al.u in contacts between a diffusive ferromagnet and a superconductor. 
Whereas in their model spin accumulation originates from a different spin conductivity of spin up 
and spin down electrons, spin accumulation originates in our model from a different interfacial 
barrier of spin up and spin down electrons. 
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